Sensitivity of Quantum Motion to Perturbation in Triangle Map 
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We study quantum Loschmidt echo, or fidelity, in the triangle map whose classical counterpart 
has linear instability and weak chaos. Numerically, three regimes of fidelity decay have been found 
with respect to the perturbation strength e. In the regime of weak perturbation, the fidelity decays 
as exp(— ce 2 t 7 ) with 7 ~ 1.7. In the regime of strong perturbation, the fidelity is approximately a 
function of et 2 ' 5 , which is predicted for the classical fidelity [G. Casati, et al, Phys. Rev. Lett. 94, 
114101 (2005)], and decays slower than power-law decay for long times. In an intermediate regime, 
the fidelity has approximately an exponential decay exp(— c'et). 

PACS numbers: 05.45.Mt, 05.45.Ac, 05.45.Pq 



I. INTRODUCTION 

The stability of quantum motion in dynamical systems, 
measured by quantum Loschmidt echo [l(, has attracted 
much attention in recent years. The echo is the over- 
lap of the evolution of the same initial state under two 
Hamiltonians with slight difference in the classical limit, 
M(t) = \m(t)\ 2 , where 

m(t) = (^ \exp(iHt/h)exp(-iH t/h)\y ) (1) 

is the fidelity amplitude. Here H and H are the unper- 
turbed and perturbed Hamiltonians, respectively, H = 
Hq + eH±, with e a small quantity and Hi a perturba- 
tion. This quantity MJt) is called fidelity in the field of 
quantum information [2J. 

Fidelity decay in quantum systems whose classical 
counterparts have strong chaos with exponential instabil- 
ity, has been studied well [3, U, la, la LD, la |9|, llfj, [llj, LL2, |l3|, 
EL EE EE E3] ■ Related to the perturbation strength, pre- 
vious investigations show the existence of at least three 
regimes of fidelity decay: (i) fn the perturbative regime 
in which the typical transition matrix element is smaller 
than the mean level spacing, the fidelity has a Gaussian 
decay, (ii) Above the perturbative regime, the fidelity has 
an exponential decay with a rate proportional to e 2 , usu- 
ally called the Fermi-golden-rule (FGR) decay of fidelity, 
(iii) Above the FGR regime is the Lyapunov regime in 
which M(t) has usually an approximate exponential de- 
cay with a perturbation-independent rate. 

Fidelity decay in regular systems with quasiperiodic 
motion in the classical limit has also attracted much at- 
tention p, [T3, [i^, liM, HM. US, HE 12^, US, H^, For 

single initial Gaussian wavepacket, the fidelity has been 
found to have initial Gaussian decay followed by power 
law decay dill, 

Meanwhile, there exists a class of system which lies be- 
tween the two classes of system mentioned above, namely, 
between chaotic systems with exponential instability and 
regular systems with quasiperiodic motion. One example 
of this class of system is the triangle map proposed by 
Casati and Prosen 26]. The map has linear instability 
with vanishing Lyapunov exponent, but can be ergodic 



and mixing with power-law decay of correlations. The 
classical Loschmidt echo in the triangle map has been 
studied recently and found behaving differently from that 
in systems with exponential instability and in systems 
with quasiperiodic motion [27J. This suggests that the 
decaying behavior of fidelity in the quantum triangle map 
may be different from that in the other two classes of sys- 
tem as well. In this paper, we present numerical results 
which confirm this expectation. 

Specifically, like in systems possessing strong chaos, in 
the triangle map three regimes of fidelity decay are found 
with respect to the perturbation strength: weak, interme- 
diate and strong. However, in each of the three regimes, 
the decaying law(s) for the fidelity in the triangle map 
has been found different from that in systems possess- 
ing strong chaos. In section II, we recall properties of 
the classical triangle map and discuss its quantization. 
Section III is devoted to numerical investigations for the 
laws of fidelity decay in the three regimes of perturbation 
strength. Conclusions are given in section IV. 



II. TRIANGLE MAP 

On the torus (r,p) € T 2 = [— 7r,7r) x [—n,n), the tri- 
angle map is 



Pn+i =p n + a sgn(r„) + 0, 

r n+l — r n + Pn+l^ 



(mod27r) 
(mod27r) (2) 



where sgn(?*) = ±1 is the sign of r for r 7^ and 
sgn(r) = for r = [26]. Rich behaviors have been 
found in the map: For rational a/ir and /3/ir, the system 
is pseudointegrable. With the choice of a = and ir- 
rational /3/ir, it is ergodic but not mixing. Interestingly, 
for incommensurate irrational values of a/V and f3/ir, the 
dynamics is ergodic and mixing. In our numerical calcu- 
lations, we take a = tt 2 and (3 = (s/E — 1)tt/2, for which 
(P/oi) is an irrational number, the golden mean divided 
by 7T, and the map is ergodic and mixing. 

The triangle map @ can be associated with the Hamil- 
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FIG. 1: (color online). Averaged fidelity at weak perturba- 
tion, a = 10 -4 (solid curve), with average taken over 50 initial 
point sources chosen randomly, N = 2 12 = 4096. The dashed- 
dotted straight line has a slope 1.7, showing that log 10 M(t) 
is approximately a function of t 1 ' 7 . For comparison, we also 
show two straight lines (dashed and dotted) with slopes 1 and 
2, respectively. 



FIG. 2: (color online). Averaged fidelity at three weak per- 
turbation strengths, a — 10 _4 (thin solid curve), 10~ 3 (dashed 
curve), and 10 -2 (thick solid curve), with average taken over 
50 initial point sources chosen randomly, N = 2 12 = 4096. 
The dashed-dotted straight line represents M 1 (t) in Eq. J5]) 
with 7 = 1.7 and c as an adjusting parameter. Inset: Fi- 
delity of a = 10~ 3 and N = 2 n ; the two curves are almost 
indistinguishable . 
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V(r) J2 S(t-nT), 



(3) 



where V(r) = — a\r\ — (3r and T is the period of kicking. 
It is easy to verify that the dynamics produced by this 
Hamiltonian gives the map @ with the replacement p = 
Tp, a — Ta, and f3 — T/3. 

The classical map can be quantized by the method 
of quantization on torus [H, H^, H3, GOJ. Schrodinger 
evolution under the Hamiltonian in Eq. ([3]) for one period 
of time is given by the Floquet operator 



U = cxp 



exp[-iV(f)], 



(4) 



where we set h = 1 in Schrodinger equation. In this quan- 
tization scheme, an effective Planck constant Tieff = T is 
introduced. It has the following relation to the dimension 
N of the Hilbert space, 



NK 



(5) 



hence, fi. e ff = 2-k/N. In what follows, for brevity, we 
will omit the subscript eff of h c g. Eigenstates of f 
and p are discretized, f\j) = jTi\j) and p\k) — kh\k), 
with j,k = -N/2,-N/2 + l,...,0,l,...,(iV/2) - 1. 
Then, making use of the above discussed relations among 
p, p, T, 5, a, /3, f3, in particular, T = h, the Floquet oper- 
ator in Eq. 0$ can be written as 



U = cxp 



exp 



~(a\r\+pr) 



(6) 



In numerical computation, the time evolution = 
U t \tpo) is calculated by the fast Fourier transform (FFT) 
method. 

The fidelity in Eq. ([1]) involves two slightly different 
Hamiltonians, unperturbed and perturbed. In this paper, 
for an unperturbed system with parameters a and /?, the 
perturbed system is given by 



+ e /?->/?. 



(7) 



Without the loss of generality, we assume e > 0. The 
parameter a — (e/h) can be used to characterize the 
strength of quantum perturbation. 



III. THREE REGIMES OF FIDELITY DECAY 

A. Weak perturbation regime 

Let us first discuss weak perturbation. As mentioned 
in the introduction, in systems with strong chaos in the 
classical limit, the fidelity has a Gaussian decay under 
sufficiently weak perturbation. The Gaussian decay is 
derived by making use of the first order perturbation 
theory for eigensolutions of H and Hq and the random 
matrix theory for AE n = E n — E® , where E n and E® are 
eigenenergies of H and Hq, respectively. Numerical re- 
sults in Ref. [HJ show agreement of the spectral statistics 
in the triangle map with the prediction of random ma- 
trix theory, hence, at first sight, Gaussian decay might be 
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FIG. 3: (color online). Variation of the averaged fidelity with 
at for a = 0. Of, 0.02 and O.f, with average taken over 100 
initial point sources chosen randomly, N = 4096. The solid 
straight line is drawn for a comparison with linear dependence 
on at. For a = 0.02 and 0.1, log 10 M(t) is approximately a 
linear function of at, before it becomes close to the saturation 
value. Inset: The distribution P(y) for the action difference 
AS at t = 40, where y = (AS - (AS))/e and (AS) is the 
average value of AS. It is calculated by taking randomly 
10 7 initial points in the phase space. P(y) does not have a 
Gaussian shape. 

expected for the fidelity decay in the weak perturbation 
regime of the triangle map. 

However, our numerical results show a non-Gaussian 
decay of fidelity for small perturbation. An example is 
given in Fig.[T]for a = i0~ 4 . To obtain relatively smooth 
curves for fidelity, average has been taken over 50 initial 
point sources (eigenstates of f) chosen randomly. This 
figure, plotted with log 10 (— log 10 M(t)) versus log 10 t, 
shows clearly that log 10 is approximately propor- 

tional to t 1 - 7 (the dashed-dotted straight line), while is 
far from the Gaussian case of t 2 and the exponential case 
of t represented by the dotted and dashed lines, respec- 
tively. 

Furthermore, we found that the averaged fidelity M(t) 
can be fitted well by 

Mi(t) = exp(-ccr 2 t 7 ) (8) 

with 7 ~ f.7 and c as a fitting parameter. In Fig. [21 we 
show fidelity decay for three different values of a. With 
the horizontal axis scaling with log 10 a 2 ^- 7 , the three 
curves corresponding to the three values of a are hardly 
distinguishable in their overlapping regions (except for 
long times). Note that, to show clearly the dashed-dotted 
straight line which represents M\(t) in Eq. we have 
deliberately adjusted a little the best-fitting value of c 
such that the dashed-dotted line is a little above the 
curves of the fidelity. 

In the inset of Fig. [21 we show curves of fidelity for the 
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FIG. 4: (color online). Fidelity decay for a = 0.1 and N = 
2", averaged over 100 initial point sources. 



same a but different values of e and N . The two curves 
are very close, supporting the assumption that e and N 
appear in the form of the single variable a as written on 
the right hand side of Eq. {8]| . This dependence of M (t) 
on the variable a for sufficiently small a can be under- 
stood in a first-order perturbation treatment of fidelity, 
as shown in the following arguments. 

Let us consider a Hilbert space with sufficiently large 
dimension N and make use of arguments similar to those 
used in Ref. H for deriving the Gaussian decay, but with- 
out assuming the applicability of the random matrix the- 
ory. It follows that, for times not very long, the averaged 
fidelity (averaged over initial states) is mainly determined 
by (exp(— iAw n t)), where Aoj n = uj n ~ uj^ and (. . .) in- 
dicates average over the quasi-spectrum. Here is an 
eigen-frequency of the Floquet operator U in Eq. ([6]) and 
uj n is the corresponding eigen-frequency of (C/e l<T ' r '). For 
large N, (exp(—iAuj n t)} can be calculated by making use 
of the distribution of Auj n . Since the two Floquet oper- 
ators U and (Ue l,jM ) differ by e m|r| , the distribution of 
Aoj n is approximately a function of a. Then, M(t) is 
approximately a function a. 

Finally, we give some remarks on the value of 7. When 
Auj n has a Gaussian distribution, M(t) has a Gaussian 
decay with 7 = 2, as in the case of systems possessing 
strong chaos. In the triangle map, the non-Gaussian de- 
cay of fidelity discussed above implies that Au n does not 
have a Gaussian distribution. Other types of distribution 
may predict values of 7 different from 2, in particular, a 
Levy distribution would give 7 < 2 in agreement with our 
numerical result. We also remark that the results here 
are not in confliction with numerical results of Ref. [13] , 
in which only the statistics of w„ (not that of Auj n ) is 
found in agreement with the prediction of random ma- 
trix theory. 
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FIG. 5: (color online). Averaged fidelity of a from 0.1 to 1, 
with average taken over 1000 randomly chosen initial pointer 
sources, TV = 2 14 = 16384. For a = 0.2 and above, the 
averaged fidelity obeys a decaying law which is different from 
that in Eq. (|9]), in particular, it is not a function of (at). 



B. Intermediate perturbation strength 



FIG. 6: (color online). Averaged fidelity at strong pertur- 
bation, with average taken over 1000 randomly chosen ini- 
tial Gaussian wavepackets, N — 2 17 = 131072. z = ei 2 ' 5 /% 
with h fixed in this figure. The solid line represents a curve 
exp(— cei 2 ' 5 ), where the fitting parameter c is determined from 
comparison with the two curves of a = 2 and 4 in the small- z 
region. 



With increasing perturbation strength, exponential de- 
cay of M(t) appears (see Fig. |3]). For a from 0.02 to 0.1, 
after some initial times and before approaching its satu- 
ration value, the fidelity decays as 



M 2 (t) = exp(-aat), 



(9) 



with a as a fitting parameter. Numerically, we found that 
a 0.08. The decay rate is proportional to (at), unlike 
in the FGR decay found in systems with strong chaos, 



M FGR (t) ~ eM-2° 2 K E t), 



(10) 



where Ke is the classical action diffusion constant [8j. 
The curves of a — 0.02 and 0.1 in Fig. [3] are quite close, 
while that of a = 0.01 has some deviation from the two. 
This implies that the exp(— aat) behavior of M(t) ap- 
pears between a = 0.01 and 0.02. Note that vertical 
shifts have been made for the two curves of a = 0.02 and 
0.1 in Fig. [3] for better comparison. 

The origin of the non-FGR decay of fidelity in this 
regime of perturbation strength, may come from weak 
chaos. In fact, in another system which also possesses 
weak chaos in the classical limit, namely, the sawtooth 
map in some parameter regime, linear dependence of the 
decaying rate on a has also been observed in the interme- 
diate perturbation regime [l3l . [Tol [33j . In this regime of 
perturbation strength, the scmiclassical theory predicts 
that, in the first order classical perturbation theory, the 
averaged fidelity is given by [l3[ 



M(t) 



dASe zAS/n P(AS) 



(11) 



where A5'(p [l , ro; t) — e J Q dt'Hi [(v(t')] is the action dif- 
ference of two the classical trajectories starting at the 
same point (po,ro) in the two systems, with H\ eval- 
uated along one of the two trajectories, and P(AS) is 
the distribution of A5(po, i"o; i). In systems possessing 
strong chaos, P(AS) may have a Gaussian form, which 
implies the FGR decay for the fidelity. In the triangle 
map, P(AS) is not a Gaussian distribution as shown in 
the inset of Fig. [3l hence, the fidelity does not have the 
FGR decay with a rate proportional to a 2 . 

It is difficult to find an analytical expression for 
P(AS), hence, we can not derive Eq. © analyti- 
cally. However, a qualitative understanding of the (at)- 
dependence of M(t) can be gained, as shown in the fol- 
lowing arguments. Equation (jlip shows that the time- 
dependence of fidelity decay comes mainly from the de- 
pendence of P(AS) on time. In the case of strong chaos, 
AS* behaves like a random walk, hence, P(AS) has a 
Gaussian form with a width increasing as \fi [8J. Since 
AS* oc e, the width of P(AS) is a function of (ey/t); then, 
Eq. (fTT|) gives the FGR decay of M(t) which depends on 
(a 2 i). In the case of the triangle map, due to the linear 
instability of the map, it may happen that the width of 
P(AS) increase linearly with t in some situations when 
t is not very long. This implies that the width of P(AS) 
may be a function of the variable (et). Then, it is possible 
for M(t) to be approximately a function of (at). 

Equation (fTTj) predicts that, up to the first order clas- 
sical perturbation theory, the dependence of M(t) on e 
and h takes the single variable a — e/Ti. Numerically 
we found that this is approximately correct, as shown in 
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FIG. 7: (color online). Averaged Fidelity for strong pertur- 
bation, from top to bottom, a = 2, 4 and 10. The average 
is taken over 1000 initial Gaussian wavepackets chosen ran- 
domly and over time from t — 2 to t + 2. N = 2 17 . The time 
axis is plotted in the logarithm scale. It shows that the long 
time decay of fidelity is slower than power law decay. 



FIG. 8: (color online). The same as in Fig. [7] with a different 
scale for the horizontal axis and for the time interval 140 < 
t < 1000. For a = 4 and 10, log 10 M(f) form two lines for 
each a. The three solid lines represent log 10 Af3(t) given by 
Eq. HU), with b = 9.6, 9.3, and 8.3 from top to bottom. 



Fig. gl Specifically, for fixed a = 0.1, M{t) of N = 2 11 
and of N = 2 12 separate at about t = 15. Indeed, for 
long times t, higher order contributions in the classical 
perturbation theory may need consideration and M(t) 
may depend on e and Ti in a different way. For larger 
iV", hence smaller h, the agreement becomes better, e.g., 
M(t) of N = 2 12 is closer to N = 2 13 than to N = 2 U . 

When a goes beyond 0.1, the exponential decay of 
M(t) expressed in Eq. © disappears, in particular, the 
dependence of M(t) on a and t does not take the form 
of (at) (see Fig. [5]). Meanwhile fluctuations of M(t) be- 
comes larger and larger with increasing a for initial point 
states. For example, Fig.[5]shows that M(t) of a = 1 has 
considerable fluctuations even after averaging over 1000 
initial point sources. Taking initial Gaussian wavepack- 
ets, the fluctuations can be much suppressed. 



C. Strong perturbation regime 

The triangle map has vanishing Lyapunov expo- 
nent, hence, its fidelity may not have the perturbation- 
independent decay which has been observed at strong 
perturbation in systems posses sing exponential instabil- 
ity in the classical limit 0, 0, Il2l Il6l |. To understand 
fidelity decay in the triangle map, it is helpful to re- 
call results about the classical fidelity given in [27]]. In 
the classical triangle map, the classical fidelity decays 
as M c i(t) ~ exp(— cet 25 ) for initial times when M c i(t) 
remains close to one, and has an exponential decay 
exp(— c'e 2 / b t) for longer times. The interesting feature 
is that the classical fidelity depends on the same scaling 
variable r = ei 2 ' 5 in different time regions. 



In the weak and intermediate perturbation regimes dis- 
cussed in the previous sections, the dependence of fidelity 
on e and t does not take the form of the single variable r. 
This is not strange, because the classical limit is achieved 
in the limit Ti — > 0, which implies a — ► oo for whatever 
small but fixed e. Therefore, it is the strong perturba- 
tion regime in which the decaying behavior of fidelity may 
have some relevance to the classical fidelity. Numerical 
results presented below indeed support this expectation. 

Figure [6] shows variation of the averaged fidelity with 
log 10 ei 2,5 , with average taken over 1000 initial Gaussian 
wavepackets chosen randomly. The initial decay of the 
fidelity of cr = 2 and 4 are quite close to the classical 
prediction exp(— eei 2 5 ). For longer times, the fidelity of 
a from 2 to 10 (with h fixed) is approximately a function 
of r, the scaling variable predicted in the classical case, 
but, the decaying behavior of fidelity is not the same 
as that of the classical fidelity, i.e., not an exponential 
decay. We found that the dependence of M(t) on % does 
not take the form of r/h, i.e., M(t) is not a function of 
the single variable (r/h). 

For long times, the fidelity has large fluctuations even 
after averaging over 1000 initial Gaussian wavepackets. 
The fluctuations can be much suppressed, when a further 
average is taken for time t . Specifically, for each time 
t, we take average over M(t') for t' from t — 2 to t + 2. 
The results are given in Fig. [3 which shows that the long 
time decay of fidelity is slower than power law decay. To 
study the decaying behavior of the slower-than-power- 
law decay, we compare it with the function 

M 3 (t)=a(log 10 i)- 6 , (12) 

with a and b as fitting parameters. In the time interval 
140 < t < 1000, the averaged fidelity can be fitted by this 
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function, as shown in Fig. where we plot log 10 M(t) 
versus log 10 (log 10 t) . Further research work is needed to 
find analytical explanations for this slower-than-power- 
law decay of fidelity. 



IV. CONCLUSIONS AND DISCUSSIONS 

We present numerical results on fidelity decay in the 
triangle map with linear instability. Three regimes of fi- 
delity decay has been found with respect to the perturba- 
tion strength: weak, intermediate and strong. At weak 
perturbation, the fidelity decays like exp(— ca 2 t 17 ). In 
the intermediate regime, the fidelity has an exponential 
decay which is approximately exp(— c'at). In the regime 
of strong perturbation, the fidelity is approximately a 
function of et 2 5 and decays slower than power law decay 



for long times. 

These results show that the fidelity in the triangle map 
obeys decaying laws which are different from those in 
systems with strong chaos or with regular motion. The 
difference is closely related to the weak-chaos feature of 
the classical triangle map. In which way and to what 
extent does weak chaos influence the fidelity decay? This 
is still an open question. Indeed, common features of 
fidelity decay in systems with weak chaos, as well as their 
explanations, should be an interesting topic for future 
research work. In particular, one may note that stretch 
exponential decay of fidelity has also been observed for 
wave packets which initially reside in the border between 
chaotic and regular regions in mixed- type systems [34| . 
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